Energy gap and proximity effect in MgB 2 superconducting wires 
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Measurements of the penetration depth A(T, H) in the presence of a DC magnetic field were 
performed in MgB-z wires. In as-prepared wires A(T, H < 130 Oe) shows a strong diamagnetic 
downturn below ~ 10 K. A DC magnetic field of 130 Oe completely suppressed the downturn. The 
data are consistent with proximity coupling to a surface Mg layer left during synthesis. A theory 
for the proximity effect in the clean limit, together with an assumed distribution of the Mg layer 
thickness, qualitatively explains the field and temperature dependence of the data. Removal of the 
Mg by chemical etching results in an exponential temperature dependence for A(T) with an energy 
gap of 2A(0)/T C » 1.54 (A(0) ~ 2.61 meV), in close agreement with recent measurements on 
commercial powders and single crystals. This minimum gap is only 44% of the BCS weak coupling 
value, implying substantial anisotropy. 

PACS numbers: 74.50.+r;74.25.Nf 



Superconducting MgB 2 [Q presents, for possibly the 
first time, a combination of phonon-mediated pairing 
together with a relatively high transition temperature 
(T c w 39.4 K ) comparable to hole-doped cuprates. Ev- 
idence for a phonon mechanism has come from several 
measurements which indicate a substantial isotope effect 

[| . Tunnelling measurements have given values of the 
energy gap ratio 5 = 2A(0)/T C ranging from 1.25 to 4 
0, B, E| fj\. NMR measurements of the n B nuclear spin- 
lattice relaxation rate give 5 ss 5 [|| while photoemission 
spectroscopy gives 5 ss 3 M. A recent tunnelling mea- 



surement has shown evidence for two energy gaps 1 10 



The temperature dependence of the London penetra- 
tion depth A is a sensitive probe of the quasiparticle den- 
sity of states and thus the minimum energy gap. Some 
early data for A in commercial MgB 2 powders showed 
apparent power law behavior suggesting nodes jll], [L2| . 
R is important to make sure that no extrinsic factors ex- 
ist that may bias the interpretation of A(T). A persistent 
complication has been the presence of surface contami- 
nants remaining from the growth process, most notably 
elemental Mg. In this paper we report magnetic screen- 
ing measurements of dense MgB 2 wires grown around a 
tungsten core. The presence of a Mg layer on as-grown 
wires gives rise to a large increase in the diamagnetic re- 
sponse below 10 K. The temperature and field dependen- 
cies of the magnetic screening are consistent with prox- 
imity induced correlations. After etching, the same wires 
show exponential behavior with a gap ratio of S « 1.54, 
less than 1/2 the BCS value of 3.53 and very close to the 
value obtained from recent penetration depth measure- 
ments on both commercial powders Jl3| , |l4[ ] and single 
crystals p5| . 

Growth of the MgB 2 wires has been described in de- 
tail elsewhere MM . In brief, boron fibers and Mg with a 



nominal ratio of MgB 2 were sealed in a Ta tube. The 
tube was sealed in quartz and placed in a box furnace at 
950°C for approximately two hours. The reaction am- 
poule was then removed from the furnace and quenched 
to room temperature. The wire samples used here had 
a tungsten core of 15 fim diameter, outer diameters of 
180 /im and 200 fim and were 2 mm long. SQUID mag- 
netometer measurements showed essentially ideal Meiss- 
ner screening (— Attx — 1) in applied fields up to 1000 Oe. 
However, tunnel diode measurements with much higher 
sensitivity revealed a clear diamagnetic downturn below 
10 K which we show was due to surface Mg. The Mg 
layer was identified by local XRD analysis and could be 
etched away with an 0.5 % solution of HCl in ethanol. 
SEM pictures after etching revealed a sinter of hexagonal 
MgB 2 crystallites with some traces of MgO. 

The penetration depth was measured with an 1 1 MHz 
tunnel-diode driven LC resonator used in several previ- 
ous studies |r7| . An external dc magnetic field (0—7 kOe) 
could be applied parallel to the ac field (~ 5 mOe) using 
a compensated superconducting solenoid. The oscillator 
frequency shift A/ = /(T) — f(T m i n ) is proportional to 
the rf susceptibility and thus to changes in the penetra- 
tion depth, AA = A(T) - A(T mi „) via A/ = -GAA, 
where G is a calibration constant |17f . The random ori- 
entation of MgB 2 crystallites implied that these mea- 
surements represent an average over in-plane and out-of- 
plane A. The polycrystalline nature also made it difficult 
to reliably estimate G for the wires. Therefore, all pene- 
tration depth data is plotted as raw frequency shift, after 
subtraction of the sample holder background. Decreasing 
frequency corresponds to increased diamagnetic screen- 
ing. 

Figure [l] shows the temperature dependence of the os- 
cillator frequency in as-prepared wires for zero DC field. 
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FIG. 1: Temperature dependence of the penetration depth at 
zero DC field. Inset: low-temperature region for two different 
diameter wires. 



The inset is a magnification of the low-temperature be- 
havior, showing a pronounced diamagnetic downturn be- 
low 10 K for two separate wires of somewhat different 
diameter. This downturn disappeared completely upon 
etching the wires. We attribute this downturn to a prox- 
imity effect induced in the Mg surface layer. Enhanced 
diamagnetism is a generic feature of proximity systems 
as carriers in the normal metal layer gradually acquire 
pairing correlations and develop a Meissner screening re- 
sponse || H p|p. 

Several quasiclassical analyses of proximity systems 
ranging from clean to intermediate have shown that the 
characteristic temperature for the appearance of screen- 
ing is 5 Ta where Ta = KVp /27rfeg<i is the Andreev 
|. For clean systems, Ta is the 



temperature 22, [23| 



temperature at which the normal metal coherence length 
£n(T) = UVp /2-Kk B T equals the normal metal layer 
thickness d. Here Vp is the Fermi velocity. For Mg the 
coherence length varies from 0.2 /im at T = 1 K to 2 /im 
at 10 K, suggesting an average Mg thickness of 2 /im 
and an Andreev temperature Ta ~ 0.6 K. 

Whether the proximity sandwich is in the clean or dirty 
limit depends upon the electronic mean free path, £ e , 
in the Mg layer, which is not known accurately. For 
example, a residual resistance ratio of 20 would give 
£ e = 0.2 /im, which must be compared to both £jv and 
d. The clean limit requires £ e >> min{^N , d} while the 
dirty limit requires £ e << £,N,d. Strictly speaking, the 
latter regime requires that both the normal metal and 
superconductor be in dirty limit, which is most likely not 
true for MgB-i p4| . Over the temperature range 1—10 K, 
all three numbers are comparable and we are likely in an 
intermediate range for which there is no analytic solution 
for the susceptibility |^4j. Uncertainties in the parame- 
ters did not justify fitting to the full numerical solutions. 



We therefore fit the data to both clean and dirty lim- 
its where analytic solutions are available in order to gain 
some qualitative understanding. 

In the clean limit the diamagnetic susceptibility of the 
normal metal layer is given by |23|] , 
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The factor 3 /4 comes from the nonlocal response in the 
normal metal layer that overscreens the external field and 
Xn(T) is a length scale given by |22|, p3|, H, 



Xn (0) 
Xn (T) 



7(A,T) 



66v (T) d/(N(T ) 
d 



(2) 



Here Xn (0) = \J Anne 2 / m ss 180 A is formally the Lon- 
don penetration depth in a superconductor with the car- 
rier mass and density of Mg. The energy gap of the 
superconductor enters through the factor, 



7 (A,T) = A/ 



irk B T+\]A 2 + {-Kk B TY 



(3) 



We take A = 2.6 meV from our own data shown later. 
For the dirty limit, a solution of the Usadel equations 
p5[ 26 leads to a power law susceptibility without any 
characteristic temperature. 
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A key feature of the proximity effect is the disappear- 
ance of screening in applied fields greater than a break- 
down field Hb (T, d) . In the clean limit this field is given 
by pi p7l, 



Hb(clean) 



V2 0o 
7T Xn (0) d 



-d./£ N (T) 



(5) 



where <po is the superconducting flux quantum and 
the result holds for T >> Ta The temperature and 
film thickness dependence of Hf,(T,d) has been verified 
in proximity systems that vary from somewhat dirty to 
clean [|o| and should therefore be applicable here. The 
breakdown field in the dirty limit is given by |E3l f27[, 



H b {dirty) w 1.9 



4>[)£e 



(6) 



Our device measures the screening of a very small AC 
field in the presence of a much larger DC field H. We 
assume that once H > Ht the AC screening vanishes. 
For H < Hi, we assume that the zero-field susceptibility 
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FIG. 2: Fits to clean limit proximity effect. Data in zero 
field were fit to obtain d ~ 1.95 fim and a ~ 1.2 am. Finite 
field (solid) curves were then generated from Eq. M. Data for 
H — 0, 20 and 130 Oe are shown. Inset: Comparison of dirty 
and clean limits at zero field. Dirty limit did not fit the data 
at any choice of parameters. 



expression holds. This approach clearly ignores nonlin- 
ear effects which a more carefully controlled experiment 
could address. The Mg layer was not uniform and we 
used a probability distribution for the normal metal layer 
thickness. The frequency shift measured upon extraction 
of the sample from the coil in-situ, combined with the ad- 
ditional diamagnetic screening, Fig. ^, gives an estimate 
of the Mg layer thickness of d « 1.62, close to the clean 
limit fit value, d 1.95. Based upon many studies of 
film growth and random processes in condensed matter 
systems we adopt a log-normal distribution of film thick- 
nesses p(x, d, a) — (y/^rrxa)^ 1 exp [(logs — d)/v2V] -2 
where d is the mean thickness and a the variance. A 
Gaussian distribution gave a less satisfactory fit. The 
proximity-enhanced diamagnetic contribution to the sig- 
nal was then taken to be, 



A/oc / X N{T)p{x:d 7 a)6{H-H b {T,x))dx (7) 
Jo 

In order to fit the data, we subtracted the signal from 
the superconductor alone, obtained after etching. (The 
superconducting signal has negligible temperature de- 
pendence below 5 K.) Data for H = were then fit 
to Eq. |?] with an overall scale factor, average thickness d 
and a as fitting parameters. (For the dirty limit we must 
also assume a mean free path.) We obtained d — 2 fim 
and a = 1.2 /j,to. These parameters were then held fixed 
and the response in a finite magnetic field was calculated. 
Fig. shows the data and generated curves for H = 0, 20, 
and 130 Oe. The finite field curves generated from the 



clean limit model all showed somewhat more screening 
than the data. This is partly due to our assumption of 
a distribution of the Mg layer thickness. Regions with 
thickness smaller than the average will have higher break- 
down fields and will continue to screen even large DC 
fields. We were not able to find a satisfactory fit to the 
dirty limit. The inset shows the best fit to the dirty 
limit that could be achieved. Given the uncertainty in 
parameters and the crudeness of the model, we feel that 
the agreement with the clean limit model is reasonably 
good. Both clean and dirty limits predict that the prox- 



N 
X 



80 



60 



40 



20 



etched 




unetched °o 10 



20 30 40 



10 



15 



20 



T(K) 



FIG. 3: Main figure: A/(T) before and after etching. Thin 
solid line shows a BCS fit as described in the text. Inset: full 
temperature scale A/(T) for etched and unetched wire. 

imity effect will exhibit substantial hysteretic effects. We 
observed no hysteresis for the range of fields shown here. 
This may be due to broadening of the first order tran- 
sition by the spread of film thicknesses. Hysteresis was 
observed at much higher fields, of order 1500 Oe. In this 
field we expect any proximity effect to be quenched but 
vortices will be present in MgB2- The hysteresis is then 
most likely due to trapped flux. We defer a discussion of 
the higher field data to another paper. 

In an effort to determine the pairing symmetry of pure 
MgB-2 the Mg layer was etched away. The result is shown 
in Fig. ||. The downturn disappeared completely and the 
temperature dependence became exponential. The inset 
shows full scale transition curves for both and etched and 
un-etched wires. The transition temperature remained 
unchanged and the only apparent change due to etching 
is the disappearance of the low-temperature diamagnetic 
downturn in A. 

Fig. U shows fits to both a full calculation (not low tem- 
perature expansion) of weak-coupling s-wave BCS form 
for A(T) and to a quadratic power law. The BCS fit gives 
a value of 2A /T C w 1.54 (2.6 meV) which is 0.43 times 
the weak-coupling BCS ratio of 3.52, implying a substan- 
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FIG. 4: Weak-coupling s-wave BCS fit (full temperature range 
calculation) and quadratic power law fit. Inset: MgBi data 
vs (T/T c ) 2 along with poly crystalline Nb for comparison. 



tial anisotropy if the material is in the weak-coupling 
limit. The inset shows the data on a (T/T c ) 2 scale along 
with data for polycrystalline Nb foil, which gives an ex- 
tremely good fit to isotropic BSC theory. The quadratic 
power law gives a poorer fit and argues against a nodal or- 
der parameter. Our value of 2.6 meV for the energy gap 
is in close agreement with recent penetration depth mea- 
surements on commercial powders |13 and single crystals 
||l5f which gave a value of 2.8 meV, and with tunnelling 
measurements which claimed a two band picture Jlpf] . 

In conclusion, we have reported measurements of the 
magnetic penetration depth in dense MgB 2 wires. We 
interpret the diamagnetic downturn in the effective A(T) 
for unetched wires as evidence for a clean limit proximity 
effect between MgB 2 and an Mg surface layer. After 
removing this Mg layer, the results are consistent with a 
minimum gap value of 2.6 meV. 
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